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Abstract. This work intends to provide a brief summary of some of tfie Galactic science issues for the next generation of very 
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about one order of magnitude more sensitivity at its central energy (at about 1 TeV), but extending the observational window 
to have a real broadband capability (from a few tens of GeV up to tens of TeV); exceeding at low energies the current VHE 
threshold for observations set by MAGIC as well as the few-tens-of-GeV sensitivity set by Fermi. Science topics regarding 
populations of emitters, pulsars and their nebula, binaries, supernova remnants, stars, and their associations, are discussed. 
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INTRODUCTION 

This work intends to provide a brief summary of some of 
the most important Galactic science issues that will focus 
our attention provided a new ground-based array, more 
sensitive, with lower energy threshold, and enhanced an- 
gular resolution as compared with the currently operating 
instruments, is available. One such forthcoming instru- 
ment may be the Cerenkov Telescope Array, or in short 
CTA, for which many European institutions have already 
joined efforts to work in a Design Study. Other papers 
in this volume (e.g., by J. Hinton, or M. Martinez) will 
give account of the current status of different aspects of 
this design. This paper is then not about on our ideas for 
CTA, nor it is a review of TeV astronomy. Because of 
that we advice the reader not to expect a thorough cov- 
erage of all current TeV results in the topics mentioned 
below. Rather, we will mention some previous results as 
a path to present ideas that could be useful in defining 
some of these future Galactic science issues; physics we 
know is there for us to test, but that still escape the reach 
of current instrumentation. 



GEV TO TEV SURVEYS: 
DISTINGUISHING POPULATIONS 

One of the first issues for Galactic science with the new 
generation of instruments will probably be to achieve a 
distinction, at a level of population studies, of the classes 
of emitters in the different energy domains. 22 VHE 
gamma-ray sources have been reported in the 2004 and 
2005 H.E.S.S. survey of the inner Galaxy (GPS). The 
third EGRET catalogue (Hartman et al. 1999, soon to be 



superseded by the advent of the new Fermi catalog!) still 
represents the companion to these VHE sources above 
an energy threshold of 100 MeV (with peak sensitivity 
between 150 and 400 MeV). It lists 271 sources, 17 of 
which are located within the H.E.S.S. GPS region. Com- 
paring the instrumental parameters of VHE instruments 
and EGRET there is a clear mismatch both in angular 
resolution and in sensitivity: In a ^ 5 hour observation 
(as a typical value in the GPS region) H.E.S.S. is a fac- 
tor of ~ 50 — 80 more sensitive (in terms of energy flux 
E^dN/dE) than EGRET above 1 GeV in the Galactic 
Plane for the exposure accumulated in the third EGRET 
catalogue. Assuming a similar energy flux output in the 
two different bands this mismatch implies at first sight 
that H.E.S.S. sources are not likely to be visible in the 
EGRET data set. Conversely (again under the assump- 
tion of equal energy flux output), VHE gamma-ray in- 
struments should be able to detect the majority of the 
EGRET sources, as has been suggested in the past. 

Figure[T]compares the energy fluxes vFv for EGRET 
sources and H.E.S.S. sources in the inner Galaxy. 
Clearly, the EGRET sources do not reach down as low in 
energy flux as the H.E.S.S. sources, a picture that should 
change with Fermi (referred to as GLAST in the Figure). 
By comparing the two surveys and analyzing the match- 
ing of the spectral properties in both energy regimes, it 
was found that (Funk et al. 2008) 

• There are rather few spatially coincident GeV- TeV 
sources for the considered Galactic region, and they 
could occur by chance, the chance probability of de- 
tecting two coincident sources within the H.E.S.S. 
GPS region is ~ 40%, thus no strong hint for a com- 
mon GeV/TeV source population is detected. 




log,o(Energy Flux (ergs cm s '^)) 



FIGURE 1. Distribution of integrated energy flux vFv for 
sources in the Inner Galaxy. For EGRET, the energy flux be- 
tween 1 GeV and 10 GeV, for the H.E.S. S. sources, the energy 
flux between 1 TeV and 10 TeV, are shown. Also shown is the 
sensitivity prediction for the GLAST-LAT for a typical location 
in the Inner Galaxy (1=10, b=0). From Funk et al. (2008). 



• Spectral compatibility (based on a power-law ex- 
trapolation) seems present for most of the position- 
ally coincident sources, but again, this is expected 
to occur by chance 

• H.E.S.S. limits at the position of the EGRET 
sources are constraining for a power-law extrapo- 
lation from the GeV to the TeV range for several 
of the EGRET sources, strongly suggesting cutoffs 
in the energy spectra of these EGRET sources in 
the unexplored region below 100 GeV. Power-law 
extrapolation of EGRET spectra seem to be ruled 
out for most of the EGRET sources. On the con- 
trary, EGRET limits at the position of the H.E.S.S. 
sources are not constraining for a power-law extrap- 
olation from the TeV to the GeV range. This picture 
may change for Fermi, which upper limits, if so are 
imposed, will be severely restricting the SEDs. 

But in addition of these facts, even from theory along 
the expectation of equal energy flux in the GeV and TeV 
band (leading to the same populations in both bands) is 
naive and can easily be wrong: EGRET sources may not 
emit comparable energy fluxes in the VHE 7-ray band 
but exhibit cut-offs or breaks in the energy band between 
EGRET and H.E.S.S.. There are well motivated physical 
reasons why the population of GeV and of TeV sources 
might be distinct: acceleration limits, particle transport 
affecting the primary spectrum of cosmic rays, convec- 
tion, particle energy losses producing steepening, par- 
ticle injection limited in time, absorption by pair pro- 
duction, redistribution of photon energies by cascading 
process. Additionally, TeV instruments are typically only 
sensitive to emission on scales smaller than ^ 1°. If any 
of the EGRET sources are extended beyond 1 ° without 
significant sub-structure on smaller scales (not precluded 
given the poor angular resolution of EGRET), current 



Imaging Cherenkov instruments may not be able to de- 
tect them since these sources would completely fill the 
field of view (FoV) and be removed by typical back- 
ground subtraction methods. A broadband instrument 
covering from tens of GeV to TeV with above Fermi sen- 
sitivity in the coiTesponding band would allow to distin- 
guish the different populations present at the high energy 
sky, separating the technical from the physical reasons 
by which these may diverge. 

PULSARS 

Pulsar science cases for VHE observatories have been 
recently rekindled by the MAGIC detection of pulsed 
emission above 25 GeV coming from Crab, revealing a 
relatively high energy cut-off in the phase-averaged spec- 
trum (Aliu et al. 2008, see the specific article presented 
in this volume with details of this measurement). This 
high energy cutoff probably indicates that emission hap- 
pens far out in the magnetosphere where the influence 
of the magnetic field, producing absorption of the form 
Y+B e^, is not enough to lead to a super-exponential 
cutoff of the high energy emission. 

Certainly, one of the fundamental unanswered ques- 
tions about pulsars concerns the mechanism and location 
by which GeV (and TeV) emission is radiated. Future 
sensitive arrays with low energy threshold could test this 
question with a large sample since essentially, pulsed de- 
tection is limited by photon statistics. Fermi, about 100 
times more sensitive than EGRET above 10 GeV, would 
allow detections more than 25 times fainter, or more than 
5 times more distant (reaching to the Galactic center). 
However, whereas Fermi sensitivity worsens with energy 
in the range £ > 10 GeV, any sensitivity of future arrays 
will significantly improve. To observationally study the 
pulsar cutoffs with a large sample, and to make a sig- 
nificant improvement in population studies (e.g., what 
happens as the observed gamma-ray pulsar luminosity 
-seemingly proportional to £'<'.65_ approaches the total 
available spin-down luminosity El how much does the 
assumption of a 1 sr beaming solid angle distort this cor- 
relation? how does this correlation look like above 10 
GeV?) is probably a mission for future ground-based ar- 
rays only. 

Especially at high energies, many Fermi sources will 
be unidentified and future aiTays could search, blindly, 
for the appearance of pulsar periods. We recall here 
that for pulsars, the phase changes as (j) = 2n{tv +t-^v), 
with t being the observation time and v = 1/P, with P 
the period and P the period derivative. The variation of 
phase due to change in the period can be neglected if 
its phase contribution is at most, less than half cycle, 
i.e., t^v < 1/2. Thus, the duration of the time for the 
search in absence of a priori information about P and P 



should satisfy Tohs, < (1 /2v)^'^. During this time, enough 
photons should be collected to determine P. To give an 
example, the useful blind observation time for Crab is 
10 hours. Thus, Fermi can only deal with the brightest 
pulsars without a-priori determined timing solutions. In- 
stead, Fermi high energy sources should all be bright de- 
tections for future VHE arrays. WiU this provide a way 
of measuring the fraction of radio quiet gamma-ray pul- 
sares directly? 

Future ground-base arrays can help complete the pul- 
sar panorama by studying extreme cases: pulsars with 
very high and very low magnetic fields, i.e., millisecond 
pulsars and magnetars. In polar cap models, the high en- 
ergy cutoff of millisecond pulsars is a decreasing func- 
tion of the magnetic field, so that pulsars with 10^ to 
10^ G could present cutoffs around or above 100 GeV, 
well in the range of future arrays, and their observations 
could constitute a direct test of the polar cap theory. On 
the other hand, the high energy emission from magnetars 
is also probably beyond the Fermi capabilities, as exem- 
plified directly in Fig. 8 of Zhang & Cheng (2002, also 
2001). Soft gamma-ray repeaters and anomalous X-ray 
pulsars are magnetars, presenting periods P = 6 — 12 s, 
being persistent X-ray sources (much beyond spin down) 
with luminosities of the order 10^"* - 10^^ erg s~\ and 
dipole fields lO'"* - 10^^ G in neutron stars. AXPs would 
be the first class of astrophysical objects with emission 
driven by magnetic field decay (Thomson & Duncan 
1996). For magnetars, we do not expect polar cap emis- 
sion: due to the large magnetic field all high energy pho- 
tons would be absorbed, i.e., in these cases, future arrays 
would be testing outer gaps models. 

Finally, what about GeV-TeV emission related with 
short-timescale pulsar phenomena? Buried beyond the 
current reach of working instruments, future arrays could 
explore any possible high-energy phenomenology re- 
lated with timing noise, in which the pulse phase and/or 
frequency of radio pulses drift stochastically, or with an- 
other type of irregularity in the radio called glitches, a 
sudden increase in the pulse frequency produced by ap- 
parent changes in the momentum of inertia of neutron 
stars. 



SNRS AND COSMIC RAYS 

Brief general considerations on current 
SNR observations: what' next? 

Measurements from current lACTs has shown that a 
few SNRs are particles accelerators up to energies above 
a few hundreds of TeV. We will exemplify with one case: 
for RX J1713.7-3946 (Aharonian et al. 2006a), photons 
up to energies £ > 40 TeV have been detected. This pro- 



vides proof of particle acceleration in the shell of RX 
J1713. 7-3946 beyond lO'^ eV, up to energies close to 
the knee in the cosmic ray spectrum. Indeed, if VHE 
gamma rays are produced via pion decay following in- 
elastic proton-proton interactions, gamma-ray energies 
of ~30 TeV imply that primary protons are accelerated 
to ~ 200 TeV in the shell of RX J1713.7-3946. On the 
other hand, if the gamma rays are due to Inverse Comp- 
ton scattering of VHE electrons, accelerated in the shell, 
off Cosmic Microwave Background photons (neglecting 
the presumably small contributions from starlight and in- 
frared photons), the electron energies at the current epoch 
can be estimated in the Thompson regime as ^ 1 10 TeV. 
Even considering Klein-Nishina effects £e ~ 100 TeV 
would be a realistic estimate. But, is a final proof for 
the gamma-ray origin (mostly leptonic or hadronic) in 
this and other sources only to come with future sensi- 
tivity, broadband coverage arrays? It must be recalled 
that Fermi, which would need about 5 years of inte- 
gration (with optimistic treatment of the diffuse GeV 
background) for a sensible distinction between the sim- 
plest hadronic and leptonic models above 1 GeV for RX 
J 17 13.7-3946 (which is in addition particularly difficult 
for Fermi since there is an EGRET source 5 times more 
luminous nearby. Funk et al. 2008), would not have an- 
gular resolution or sensitivity enough to do this for most 
other (all others?) SNRs, and it will certainly be very dif- 
ficult for it to determine space dependent changes in the 
spectrum of the emitted gamma-rays in the way it has 
been already determined by current instruments. 



Observing non-linearity in shocks 

Fermi 1st order acceleration is probably not a good 
enough description of particle behavior in SNR shocks. 
In Fermi theory, typical power law spectra with slope 
~ 2 — 2.4 are generated as a result of the acceleration 
process. But even when a small fraction of particles 
are injected in the acceleration box, they can take away 
an appreciable fraction of the available energy, making 
the dynamical reaction of the accelerated particles non- 
negligible (for a review see MaUcov & Drury 2001). If 
so, the theory of diffusive particle acceleration at shocks 
becomes non-linear, predicting a high acceleration effi- 
ciency, a magnetic field amplification, and particle spec- 
tra that is no longer described simply by power-laws. Fu- 
ture broadband observations of SNRs could find these 
effects. 



Diffusion 

Perhaps one of the secure science cases for cosmic- 
ray physics with future arrays involve the study of dif- 
fusion. While travelling from the accelerator to the tar- 
get the spectrum of cosmic rays is a strong function 
of time f, distance to the source R, and the (energy- 
dependent) diffusion coefficient D{E) (e.g., Aharonian & 
Atoyan 2006). Depending on these parameters t,R,D{E) 
one may expect any proton, and therefore gamma-ray 
spectrum: hard/soft/with and without TeV tail/with and 
without GeV counterpart, etc., and it gets even worse if 
the target is moving, like in the case of a stellar wind, 
where the process of convection can also play a role. Fu- 
ture instruments will open the parameter space in these 
three quantities, which is mostly unaccessible to current 
experiments. A determination, with high sensitivity, of 
spatially differentiated gamma-ray sources related to the 
same accelerator would lead to the experimental deter- 
mination of the local diffusion coefficient D{E) and/or 
of the local injection spectrum of cosmic rays (see, e.g., 
Gabici and Aharonian 2007). 

Perhaps one of the best candidates for a diffusion 
of cosmic-rays scenario observed already with the cur- 
rent generation of instruments is IC443. Discovered by 
MAGIC (Albert et al. 2007b), the steep-spectrum VHE 
source is spatially dislocated from a nearby GeV detec- 
tion. The distribution of molecular clouds and the multi- 
frequency phenomenology are consistent with the inter- 
pretation of cosmic-ray interactions with a giant molec- 
ular cloud lying in front of the remnant, producing VHE 
gamma-rays but no counterpart at lower energies. At 
a morphological level, then, the lower the energy, the 
more coincident with the SNR the radiation should be 
detected. At a spectral level, sufficient statistics should 
show that the lower the gamma-ray energy the harder the 
spectrum is (it is ~ —2 in the EGRET range). Both pre- 
dictions should show in future observations, and with as 
a combination of GeV and TeV data. Figure |2] provides 
a summary of this case. It is a clear tip of the iceberg for 
cases using future arrays, which broadband expected ca- 
pability, and enhanced sensitivity will provide an arena 
for immediate testing of these ideas. 

In the framework of diffusion models, Figure[3]shows, 
as contour plots, the energy at which the maximum of 
the SED is to be found for the cases of an impulsive ac- 
celeration of cosmic rays, at different distances, ages of 
the accelerator, and diffusion coefficient. The diffusion 
radius, for times less than the proton-proton interaction 
timescale, t <^ Tpp, is RdifiE) — 2y^D{E)t, so that at a 
fixed age and distance, only particles of higher energy 
will be able to compensate a smaller Z)io, producing SED 
maxima at higher /i-values. Small values of Dio are ex- 
pected in dense regions of ISM. It is interesting to note 
that for many, albeit not for all, of the SEDs, the maxi- 



mum in E^F space is found at energies beyond the en- 
ergetic range of Fermi, opening thus the gate for stud- 
ies with the next generation of instruments, particularly 
observing in the few tens of GeV. An observational dis- 
covery of a 1 - 100 GeV maximum, provides interest- 
ing clues about the nature of the astrophysical system 
that generates the gamma-rays: the range of accelerator- 
target separations and ages of accelerator that would pro- 
duce such a 1-100 GeV maximum is rather limited (see 
in Figure|3]the narrow contours for maxima at such ener- 
gies), which would help to produce a direct identification 
of the source, in case such system is found in the vicinity 
of the detection. 



Molecular clouds 

Another related issue concerns the study of cosmic-ray 
diffusive penetration into molecular clouds. If the dif- 
fusion coefficient inside a cloud is significantly smaller 
as compared to the average one in the neighborhood, 
low energy cosmic-rays are excluded from penetrating 
deep in the cloud, and part of the gamma-ray emission 
from the cloud is suppressed, with the consequence of its 
gamma ray spectrum appearing harder than the cosmic- 
ray spectrum (see e.g., Gabici et al. 2007). Both of this 
effects are to be more pronounced in the densest central 
region of the cloud. Thus, an angular resolution of the or- 
der of 1 arcmin or less would allow to resolve the inner 
part of the clouds and measure the degree of penetration 
of cosmic-rays into them, allowing to put constraints on 
the change of the diffusion coefficient. 

In addition, future instruments will be sensitive 
enough as to detect even with no cosmic-ray enhance- 
ment, a population of close (100-300 pc) small clouds, 
most of them translucent^ These were discovered in 
recent sensitive CO surveys at relatively high Galactic 
latitudes, e.g., covering \b\ < 30 deg and 5 > -17 deg 
(1 < 230 deg) with a sampling interval of 1/4 deg or 
better (Dame & Thaddeus 2004). There are about 200 
relatively small and isolated molecular clouds \b\ > 10 
degrees, with masses between 1 and 100 solar. Many of 
these will be detectable by future-ACTs, what will allow 
us to feedback upon the mass estimation of each of the 



Small molecular clouds are divided into three classes on the basis 
of their visual extinction and astrochemical properties (van Dishoeck 
& Black 1988): diffuse, translucent, and dark. Translucent clouds, 
the intermediate regime, have CO abundances in the range 10^^ — 
lO^"* and column densities greater than lO" cm^^. Different from 
dai'k clouds, where the chemistry is driven by collisional processes, 
translucent clouds are still dominated by photoprocesses. It is in this 
translucent regime where most of the carbon becomes molecular, and 
so they are detectable in CO, but are relatively optically thin - and thus 
their name - exhibiting low optical extinctions. 
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FIGURE 2. Left: Sky map of y-ray candidate events (background subtracted) in the direction of MAGIC J0616+225 for an 
energy threshold of about 150 GeV in galactic coordinates. Overlayed are '^CO emission contours (cyan), contours of 20 cm VLA 
radio data from (green), X-ray contours from Rosat (purple) and y-ray contours from EGRET (black). The '^CO contours are at 7 
and 14 K km/s, integrated from -20 to 20 km/s in velocity, the range that best defines the molecular cloud associated with IC 443. 
The contours of the radio emission are at 5 mjy/beam, chosen for best showing both the SNR IC 443. The X-ray contours are at 
700 and 1200 counts / 6 ■ 10^^ sr. The EGRET contours represent a 68% and 95% statistical probability that a single source lies 
within the given contour. The white star denotes the position of the pulsar CXOU J061705. 3-1-222127. The black dot shows the 
position of the 1720 MHz OH maser . The white circle shows the MAGIC PSF of tr = 0.1°. From Albert et al. (2007b), see also 
that paper for references to the multi-frequency data. Right: MAGIC and EGRET measurement of the neighborhood of IC 443 
(stars and squares, respectively) as compared with predictions of a model where cosmic-ray diffuse away from the 30 kyr SNR IC 
443. At the MAGIC energy range, the curves show the predictions for a cloud of 8000 Mq located at 20 (1), 25 (2), and 30 (3) 
pc, whereas they correspond to 15 (1), 20 (2), 25 (3), and 30 (4) pc in the bottom panel. At the EGRET energy range, the curve 
shows the prediction for a few hundred Mq located at 3^ pc. The EGRET sensitivity curve (in red) is shown for the whole lifetime 
of the mission for the Galactic anti-centre (solid), which received the largest exposure time and has a lower level of diffuse y-ray 
emission, and for a typical position in the Inner Galaxy, dominated by diffuse y-ray background. The Femi sensitivity curve (in blue) 
(taken from http://www-glast.slac.stanford.edu/software/IS/glast_latperformance.htmll show the 1-year sky-survey sensitivity for 
the same position in the Inner Galaxy. From Torres et al. (2008). 
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FIGURE 3. For each combination of age and accelerator-target separation, for which more than two thousand spectra where 
numerically produced, the energy of the maximum of such spectra are shown in a contour plot. The color of the different contours 
corresponds to the range of energy where the maximum is found according to the color bar on the top of the figures. From left 
to right, plots are created for the case of an impulsive source injecting protons in a medium with Dio=10^^cm^/s, lO^^cm^/s and 
lO^^cm^/s. After Rodriguez et al. (2008). 



clouds, and/or the variation of cosmic-ray spectrum in 
the Earth vicinity (Torres et al. 2005). 



Galactic Center 

Having said all above, it is straightforward that the 
Galactic Center itself (as discussed in more detail in 
other contributions to this volume) will be one of the 
prime science targets for the next generation of VHE in- 
struments. Essentially, following the results by Aharo- 
nian et al. 2004, 2006b; Albert et al. 2006b), one can 
realize of the need of better sensitivity and angular res- 
olution to be able to pursue an analysis of cosmic-ray 
diffusion in the region, with a few pc binning (implying 
about 1 arcmin of angular separation), distinguishing be- 
tween one or multiple originators of the primary cosmic- 
ray population. 

Indeed, apart from the candidate counterparts found in 
the error circles of the measured detections of the central 
Galactic Center source observed both by MAGIC and 
HESS: the shell-type SNR Sgr A East, the pulsar wind 
nebula G 359.95—0.04, and the supermassive black hole 
Sgr A* itself; as well as the nearby TeV source coincident 
with the SNR G 0.9+0. 1, it is interesting to focus on what 
lies beneath. The diffuse Galactic center emission turned 
out to be visible after subtraction of the previous two TeV 
sources, as reported by Aharonian et al. (2006b). These 
diffuse emission presented two main features: it appears 
correlated with molecular clouds (as for instance, traced 
by the CS (1-0) hne), and it by far exceeds (by a fac- 
tor of 3 to 9) the gamma-ray emission that would be 
produced if the same target material would be subject 
to the local neighborhood cosmic-ray environment. The 
observed morphology and spectrum of the gamma-ray 
emission provide evidence that one or more cosmic-ray 
accelerators have been active in the Galactic center in the 
last 10000 years. These accelerator, then, could provide 
the cosmic -rays that after hadronic interactions produce 
the diffuse emission seen. An alternative (which future 
instruments should rule out (or not), is to entertain the 
possibility for a number of Inverse Compton sources to 
be behind the emission (albeit the reasons for a correla- 
tion with the molecular clouds would not be that clear in 
such a case). 

The ratio of gamma-ray emission to molecular density 
varies with galactic longitude in a uniform way, except 
for a pronounced dip at /« 1.3° in the TeV emission 
that is not mimicked by the molecular content (Aharo- 
nian et al. 2006b). The implied non-uniform cosmic ray 
density can be explained if cosmic rays injected close to 
Sgr A have not yet had time to diffuse out to the / w 1 .3° 
region. Assuming a fixed value for the diffusion coeffi- 
cient of £)= 10^° cm s , what is typical for the Galac- 



tic Disk at TeV energies, an angular displacement of 1 ° at 
the distance of the Galactic Center corresponds to a dif- 
fusion time of ~10"* years, close to the age of the SNR 
Sgr A East. Future instruments with enhanced sensitiv- 
ity, angular resolution, and energy coverage would allow 
to go deeper into this cosmic-ray propagation issue, and 
actually determine the energy dependent diffusion coef- 
ficient observationally (e.g., Hinton et al. 2006). A larger 
sensitivity would also allow similar kind of studies to be 
done in other Galactic regions, closer to Earth, where the 
cosmic ray environment is probably less enhanced, and 
that are also beyond the reach of current lACTs. 



PWN 

As it is well known, the Crab Nebula is a very ef- 
fective accelerator (exemplified from observing radiation 
across more than 15 decades of energy) but not an effec- 
tive inverse Compton gamma-ray enutter. Indeed, we see 
gamma-rays from Crab because of its large spin-down 
reservoir (^ 10^^ erg s^'), although the gamma-ray lu- 
minosity is much less than its spin-down power, what 
can be understood because of a large magnetic field, 
whose strength also depends on the spin-down reservoir. 
A less powerful pulsar would imply a weaker magnetic 
field, what would imply a higher gamma-ray efficiency 
(i.e., a more efficient sharing between synchrotron and 
inverse Compton losses). For instance, HESSJ 1825- 137 
(Aharonian el al 2006c) presents a similar luminosity to 
Crab, but has 2 order of magnitude less spin down power, 
and its magnetic field has been constrained to be in the 
range of a few /iG, instead of hundreds. The differen- 
tial gamma-ray spectrum of the whole emission region 
from this latter object is measured over more than two or- 
ders of magnitude, from 270 GeV to 35 TeV, and shows 
indications for a deviation from a pure power law that 
future arrays could observe. Forecasting future obser- 
vations, spectra have also been determined for spatially 
separated regions of HESS J1825-137. The photon in- 
dices from a power-law fit in the different regions show a 
softening of the spectrum with increasing distance from 
the pulsar and therefore an energy dependent morphol- 
ogy. If the emission is inverse Compton generated, the 
pulsar power is not enough to generate the gamma-ray 
luminosity, suggesting that the pulsar had a higher injec- 
tion power in past. Is this common for other PWN? What 
can that teach us about the evolution of pulsar winds? 
Another notable case is that of Vela X (Aharonian el al 
2006d), for which the first detection of what appears to 
be a VHE Compton peak in the SED diagram was found. 
Although a hadronic interpretation has also been put for- 
ward (Horns et al. 2006) it is yet unclear how large the 
content of ions in the pulsar wind could be. 



All in all, the study of pulsar wind nebulae (PWN) 
with future arrays opens the possibility to focus on sev- 
eral issues related to pulsar winds, and how they reflect 
on their capability of cosmic-ray accelerators: magneti- 
zation, composition, bulk Lorentz factor, anisotropy, par- 
ticle distributior0. 



Absorption by the local ISRF 

Given the known existence of multi-TeV sources, fu- 
ture arrays could also measure VHE absorption in the 
interstellar radiation field (ISRF). This is in fact impossi- 
ble for other experiments, like Fermi, due to insufficient 
energy coverage, and very hard or impossible for current 
experiments, due to insufficient sensitivity. The attenua- 
tion due to the CMB is at only 10 Kpc if E >500 TeV. 
But the attenuation due to the ISRF (which has a com- 
parable number density at longitudes 20/xm to 300/xm, 
can produce absorption at about 50 TeV (Zhang et al. 
2006, Moskalenko et al. 2006). Observation on the cut- 
off energy of different sources will provide independent 
information to test and constrain the ISRF model. 



MICROQUASARS AND GAMMA-RAY 
BINARIES 

Very recently, a few massive binaries have been iden- 
tified as variable VHE 7-ray sources. They are PSR 
B1259-63 (Ahai-onian et al. 2005a), LS 5039 (Aharonian 
et al. 2005b, 2006), LS I H-61 303 (Albert et al. 2006, 
2008a,b), and Cyg X-1 (Albert et al. 2007). The nature 
of only two of these binaries is considered known: PSR 
B 1259-63 is formed with a pulsar whereas Cyg X-1 is 
formed with a black hole compact object. The nature of 
the two remaining systems is under discussion. 

The high-energy phenomenology of Cyg X-1 is dif- 
ferent from that of the others. It has been detected just 
once in a flare state for which a duty cycle is yet un- 
known. The three other sources, instead, present a be- 
havior that is fully correlated with the orbital period. The 
latter varies from about 4 days in the case of LS 5039 to 
several years in the case of PSR B 1259-63. Cyg X-1 and 
the the remaining TeV sources also differ in their SEDs. 
In Cyg X-1, the transient VHE luminosity was less than 
1% of the X-ray luminosity. For PSR B 1259-63, LS 5039 
and LSI H-61 303, the emission above 10 MeV dominates 
the radiative output. Are these differences indicative of a 
universal behavior related with the system composition? 
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See more about this in tlie section about binaries. 



FIGURE 4. Lightcurve for LS 5039, together with power- 
law fittings in 0.1 -phase intervals. From Aharonian et al. 
(2006). 



Figure IH from Aharonian et al. (2006), where the 
reader can also check the spectral variability, presents the 
current public data on LS 5039. It is one of the richest 
datasets of high energy astrophysics, and certainly in the 
top tier of Galactic ones. In the top of the Figure, the y- 
ray flux (F > \ TeV) lightcurve (phasogram) of LS 5039 
from HESS data (2004 to 2005) is presented on on a run- 
by-run basis. Each run is '--^28 minutes. Two full phase 
periods are shown for clarity. The blue solid arrows de- 
pict periastron and apastron. The thin red dashed lines 
represent the superior and inferior conjunctions of the 
compact object, and the thick red line depicts the Lomb- 
Scargle Sine coefficients for the period giving the high- 
est Lomb-Scargle power in the periodicity analysis, re- 
sulting in the TeV determination of the orbital period. 
These data has been binned in 0.1 -phase intervals, and 
subsequently fitted with a power law, which normaliza- 
tion and index are shown in the middle and bottom panel 
of the figure (for energies 0.2 to 5 TeV). Because of low 
statistics, and presumably too, of vastly different statis- 
tics at the higher end of the spectrum in each of these 
bins, more complicated functions such as a power-law 
with exponential cutoff provide a no better than a pure 



power-law. This is then a concept for showing a power- 
ful testing machinery: the more sensitive the instrument, 
the smaller the phase bin (the interval) upon which an 
spectrum can be given. Detailed models should then aim 
in providing spectral evolution along this and others bi- 
nary orbits, taking into account the most complete possi- 
ble physical (e.g., geometrical effects, anisotropic fields, 
cascading). 

The referred question about what kind of systems are 
these, now remains open. Even, it might be possible 
to say that it has been worsened since the Swift-BAT 
detection of a SGR/AXP-like burst from the direction 
of LS lH-61 303 (see the GCN circular 8209 and their 
foUow ups). The ascription of this short 0.2 s-burst to LS 
1 +61 303 itself, what would qualify it most likely as a 
magnetar (see e.g., ATEL 1715) has been put to question 
by several authors (see ATELs 1730, 173 1, 1740). Recent 
reports of observations made by RXTE have also called 
to surprise, finding that the X-ray emission from LS 1 H-6 1 
303 is changing by up to a factor of six over timescales of 
several hundred seconds as well as doubUng in times as 
fast as 2 seconds (Smith et al. 2008). We recall, however, 
that RXTE has a large position uncertainty, and since 
no follow-up observations to any of these flares were 
made, it is perfectly possible that they were generated 
in a source nearby (in sky projection) but different from 
LS I-H61 303. 

We are not sure whether the identification of sources 
such as LS 5039 or LS I +61 303 will come before or 
after the existence of a new generation of ground based 
gamma-ray instruments, despite current efforts. The ul- 
timate proof of a pulsar in these systems, if such ex- 
ist, i.e., the detection of pulsed emission at any wave- 
length, might remain obscured by strong absorption (see 
e.g., Zdiarski et al. 2008). If so, in order to confirm or 
rule out a particular composition for these or other sys- 
tems, it might be necessary to recur to high energy astro- 
physics, and compare the most detailed predictions pos- 
sible with observations having enough statistics in short 
time-intervals. 

If gamma-ray binaries are pulsars, is the gamma-ray 
emission coming mostly from processes within the pul- 
sar wind zone or rather from those being produced by 
particles accelerated in the wind collision shock? Is this 
answer a function of energy? Gamma-ray astronomy can 
indeed provide clues not only to answer this question, 
but also, to acquire information regarding the particle's 
energy distribution within the pulsar wind zone itself. 
For instance, data on LS 5039, analyzed in the context 
of pulsar wind zone models already rule out that mono- 
energetic electrons are responsible for most of the emis- 
sion (Sierpowska-Bartosik 2008, Cerruti et al. 2008). 
What else can future data tell us? How would that im- 
pact on current models for dissipation in pulsar winds 
(see, e.g., Jaroschek et al. 2008 and references therein)? 



It is not implausible that close systems may trigger dif- 
ferent phenomenology within the pulsar wind zone, ul- 
timately affecting particle acceleration there. Models of 
particle energization and dissipation in pulsar winds are 
currently made for isolated objects, it is probably true to 
say that we have no clue whether the inclusion of such 
objects in close binaries will affect the wind behavior 
or even the magnetosphere in appreciable ways. Could 
we gain knowledge on this using future short-timescales 
gamma-ray observations? 

Among the possibilities for future instruments, it is 
also worth noticing that of the determination of the duty 
cycles of high energy phenomena. For example, contin- 
ued observations of key objects (such as Cyg X-1) with 
current instrument's sensitivity, using sub-arrays of fu- 
ture ACTs provides a secure science case. A constant 
coverage able to measure 0.1 Crab fluxes in less than 
1-hour flares, at a Galactic center distance, can be used 
to measure the duty cycle and provide triggering condi- 
tions for the full array of the same TeV instruments, as 
well as to others at different frequencies (trigger and be 
trigger by, e.g., LOFAR). We recall that energy coverage 
in 10-100 GeV band, in the case of short flares, is not 
possible with current instruments (AGILE or Fermi lack 
sensitivity); nor it is practically possible for higher ener- 
gies, since that would imply many-months to year-long 
campaigns observing the same objects. 

Improved angular separation at high energies will also 
provide new science cases for microquasars, particu- 
larly if their jets contain a sizable fraction of relativistic 
hadrons. While inner engines will still be out of domain 
for future arrays of high energy telescopes, the study of 
any possible microquasar jet/ISM interaction might not 
be, leading to distinction between the gamma-ray emis- 
sion produced by the central object (which may be vari- 
able) and that produced by the ISM interaction (which 
may be stable, as in the case of cloud overtaken by a SNR 
shock), see e.g., Bosch-Ramon et al. (2005). To make this 
distinction possible, at typical kpc distances from Earth, 
few arcmin from the compact object should be separated. 

Finally, also in the case of microquasars, it is known 
that black holes display different X-ray spectral states: 
with transitions between a low/hard state, where a com- 
pact radio jet is found, to a high/soft state, where the 
radio emission is quenched by large factors, or there 
is no detectable radio emission at all (see e.g.. Fender 
2003). Are these spectral changes related with gamma- 
ray emission? Is there any gamma-ray emission during 
non-thermal radio flares (which can reach up to a factor 
of a 1000 increase in radio flux density)? Indeed, gamma- 
ray through inverse Compton is expected when the radio- 
to-X-ray emission during the flares is modeled via syn- 
chrotron radiation of relativistic electrons suffering ra- 
diative, adiabatic and energy-dependent escape losses in 
fast-expanding plasmoids (radio clouds) (see e.g.. Aha- 



ronian & Atoyan 1999). Can future gamma-ray observa- 
tions put constrains upon tlie magnetic field in the plas- 
moid? Needless to say, the cross-correlations of large- 
coverage instruments in the radio (e.g., LOFAR) and 
gamma-rays will be essential in tackling these issues. 

STARS AND STAR ASSOCIATIONS 

O and B stars lose a significant fraction of their mass in 
stellar winds with terminal velocities of order Voo = 10^ 
km s^'. With mass loss rates as high as — (10^'' — 
10^^) Mq yr^\ the density at the base of the wind can 
reach 10^'^ g cm^^. WR stars represent an evolved stage 
of hot (T^ff > 20000 K), massive (M^amS > ^5Mq) 
stars, and display some of the strongest sustained winds 
among galactic objects. Thus, colliding winds of mas- 
sive star binary systems involving WR, O, or B stars are 
interesting candidates to generate gamma-ray emission. 
When two of these massive stars are in binary systems, 
the region of the hydrodynamical equilibrium of their 
winds forms a shock that can accelerate particles: non- 
thermal spectra in the radio band (e.g., as in the case of 
Cyg OB2 n.5, or WR 147, e.g., see Wilhams et al. 1997) 
have been accordingly detected, as well as hints of cor- 
relation of yet unidentified EGRET gamma-ray sources 
with such binaries have been found (e.g., Romero et al. 
1999). 

With different levels of detail, leptonic (inverse Comp- 
ton of relativistic electrons with the dense photospheric 
stellar radiation fields in the wind-wind colUsion zone, 
e.g., see Reimer et al. 2006 and references therein) and 
hadronic (neutral pion decay products, where mesons 
produced by inelastic interactions of relativistic nucle- 
ons with the wind material produce the gamma-rays, e.g., 
see Benaglia et al. 2001) scenarios have been developed. 
Additionally, inverse Compton pair cascades initiated by 
high-energy neutral pion decay photons (from nucleon- 
nucleon interactions in the stellar winds, e.g., Bednarek 
2005) or collective wind scenario in young stellar cluster 
or OB-associations (e.g., diffusive shock acceleration by 
encountering multiple shocks, e.g., Klepach et al. 2000; 
or inelastic proton interactions with collective winds tak- 
ing into account the possible convection of cosmic ray 
primaries, e.g., Torres et al. 2004a, Domingo-Santamaria 
& Torres 2006) have been also put forward. The sensi- 
tivity of current lACTs is not enough to test them thor- 
oughly and thus these models are being developed with 
basically no observational feedback at the highest ener- 
gies. 

Indeed, sensitive upper bounds upon the gamma-ray 
emission from isolated binary systems (to avoid the com- 
plexity in the determination of origin of the radiation be- 
tween binaries and collective effects, see the discussion 
below) were just recently presented by the MAGIC tele- 



scope (Albert et al. 2008b). Thus the establishment of 
WR binaries as VHE y-ray sources is yet pending. The 
MAGIC observations of the system WR 147 show al- 
ready that Fermi should see a flux cutoff well within its 
range of detectability in the tens of GeV regime, if it is 
able to detect these stars at all. Within the system param- 
eters mentioned due to the low expected opacity to y-ray 
escape, MAGIC observations would directly limit the ac- 
celeration efficiency in this system so as not to produce 
a significant number of particles with energies typically 
greater than 100 GeV. The next generation of instruments 
will allow us to go from limits to detection, with the con- 
sequent possibility of testing physical assumptions. 

Since the HEGRA measurements of the Cygnus OB2 
neigborhood and the detection of TeV J2032+41 30 (Aha- 
ronian et al. 2005), which is still unidentified (Albert et 
al. 2008), little doubt was left that there exists a direct 
connection between high-energy photon emission and re- 
gions of star formation as a whole. The observational 
connection referred to has recently been more firmly es- 
tablished with the detection of the Wd2 cluster (Aharo- 
nian et al. 2007). This contains a population of at least 8 
stars earlier then 07, 2 WRs, and in particular WR20a, 
the most massive measured star (8 1 solar masses) in our 
Galaxy (a WN6+WN6 binary). This H.E.S.S. source, 
that shows no variability nor periodicity and is extended, 
is also still unidentified: it is not yet known whether the 
emission is produced by a single or by several cluster 
member (say WR20a and others) or by all, in collective 
effects within the cluster Unless there are extreme dif- 
ferences in the spatial extent of the particle distributions 
producing radio. X-ray, and VHE gamma-ray emission, 
scenarios based on the colliding stellar winds in the WR 
20a binary system face the problem of accounting for a 
source extension of 0.2 degrees in the VHE waveband. 
But despite this, unless an orbital period is detected in 
future datasets, what in fact will not generically be pos- 
sibly due to the fact that most expected TeV emitters are 
long, multi-year period binaries, it would be very difficult 
or impossible for the current generation of instruments to 
distinguish whether the radiation observed from these as- 
sociations (Wdl is an example, but there is no reason to 
suppose Wd2, 30 Dor in the LMC, and other associations 
will not be detected as well) is coming from a few iso- 
lated binary systems or is generated as a collective effect 
of the whole cluster. Rather, the way to make this distinc- 
tion lies perhaps in the combined power of an improved 
angular resolution and sensitivity of future instruments, 
which will allow to produce detailed maps of gamma- 
ray emission at different energy cuts in order to compare 
with theoretical predictions. 

Finally, future instruments can aim constructing an 
overaU understanding of the relationship between star 
formation processes and gamma-ray emission. They can 
experimentally establish whether there is a direct corre- 



lation between star formation rate and gamma-ray lumi- 
nosity when convection and absorption processes at the 
different environments are taken into account. If such 
correlation exist, how does it change with energy? In- 
deed, the discussed stellar associations as Galactic can- 
didates for very high energy emission immediately find 
extragalactic counterparts in starbursts and luminous in- 
frared galaxies. For the latter, we currently have only 
upper limits for the nearest starburst (NGC253, Aha- 
ronian et al. 2005c) and the nearest ultra-luminous in- 
frared galaxy (Arp 220, the place in the universe with 
the highest supernova explosion rate known, Albert et 
al. 2007c). These galaxies are expected to appear in 
gamma-rays provided enough sensitivity is achieved. De- 
tailed models have been constructed for NGC 253 which 
could be tested already with the current instruments (e.g., 
Domingo-Santamaria & Torres 2005); others, like for in- 
stance the case of Arp 220 (Torres 2004) require yet fur- 
ther sensitivity. 

CONCLUDING REMARK 

Galactic physics will blossom with a new-generation of 
instruments like the planned CTA, and after the pathmak- 
ing findings of the Fermi survey. It will enter into an era 
of precision high-energy observations, it will not be just 
a discovery trip, as it has partly been up to now, but a 
real immersion into understanding the physical processes 
generating the emission detected and into detailed model 
testing. 
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